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Abstract 

In this paper, we give a numerical method for pricing long maturity, 
path dependent options by using the Markov property for each underlying 
asset. This enables us to approximate a path dependent option by using 
some kinds of plain vanillas. We give some examples whose underlying 
assets behave as some popular Levy processes. Moreover, we give some 
payoffs and functions used to approximate them. 
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1 Introduction 

In this paper, we give a numerical method for pricing some path dependent 
options by using the Markov property for each underlying asset process. Path 
dependent options are options whose payoff at maturity depend on the past 
history of the underlying asset as well as the price at maturity. Asian options, 
Lookback options and barrier options are thier well-known examples. Some of 
them as the above are difficult to calculate analytically and numerically, and 
there have been numerous studies on how to do this. 

One numerical method of pricing path dependent options uses Monte Carlo 
simulation. It has an advantage in that with it, we can simulate the expectation 
of the payoff without detailed discussion of the payoff type, but only with the 
distribution of the underlying asset. N. Hilber, N. Reich, C. Schwab and C. 
Winter [TT] or Cont, R and Tankov, P [7] give the Monte Carlo method when 
market models are extended to Levy processes. They mention that most of 
Levy processes can be simulated approximately. As seen from them, the Monte 
Carlo method has the applicability for many payoff types and many underlying 
asset models. However, the longer a maturity time is, the more the calculation 
time is consumed. 

Several analytical methods of path dependent options have been also devel- 
oped. An analytical formula for standard barrier options is given in papers by 
Merton [T7] or Reiner and Rubinstein [TH] when the underlying asset behaves 
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as the geometric Brownian motion. An analytical formula for partial barrier 
options is also given by Armstrong [3] or Heynen and Kat [10] under the same 
conditions. When the underlying asset is more general such as a Levy pro- 
cess, it is not easy to find the exact value of them. There are some studies 
about this problem, for example, Kudryavstev and Levendorski |13j give the 
Fast Wiener-Hopf factorization method. Asian options are difficult to calculate 
analytically. Geman and Yor [9] give a semi-analytical formula by using the 
Laplace transformation. Linetsky [15] uses a spectral expansion approach, and 
Benhamou [I] uses a convolution method for pricing them. Albrecher and Pre- 
dota [5] approximate the arithmetic option price based on the moments of the 
average. 

Most approaches of these methods depend on the payoff functions. In con- 
trast, a method using the Malliavin-Watanabe calculus given in Kumitomo and 
Takahashi [14] is valid for some path dependent options. Bermin [6] and [5] show 
that the Malliavin calculus approach can be applied for any square integrable 
payoff. 

Our method is an analytical approximation approach. In [12j . we give the 
asymptotic behavior of the prices of an Asian option about maturity time T 
as the underlying asset behaves as the geometric Brownian motion. When the 
maturity time T — * oo, the expectation of the payoff for an Asian option can 
be approximated by that of the payoff for some linear combinations of plain 
vanillas; i.e., there exists a non-zero constant, D, such that when T — > oo, 

N N 

E ^N £ St ^ - *0 + l ~nJ2 E[(Sr-T t - K) + ] + Da{T) 

i=l i=l 

where the underlying asset is defined by S t = Sq exp(a\V t + (r — j(j 2 )t), If is a 
strike price, {t{\ is fixed time, W = {Wt}a<t<T is a one-dimensional Brownian 
motion, and a is defined by 

(1-1) a(T) = -Lex P (-^(,-V)r). 

In this paper, we generalize the above idea to an almost universal situation; i.e., 
to some payoff functions and a class of the underlying asset process. 

To summarize our result of this paper, let the underlying asset behave as 
positive Markov process St- We represent an expectation of a path dependent 
option with maturity time T as 

A*(T)=E[g((S s ) T - T < s<T )) 

where g is a function from path space which represents the payoff of the path 
dependent option. Roughly speaking, if St decays uniformly in M; i.e., there 
exists a* : [0, oo) — + (0, oo) such that P ^f^ "convergence" to non-trivial mea- 
sure as T — > oo, and also if g satisfys "good" conditions, there exists A*(T) and 
a constant, C, such that the error term between A* and A* can be estimated 
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when T — > oo; i.e., 

^(T)-1*(T) 

In practical sense, if A*(T) is either easy to calculate or is quoted in the market, 
we easily obtain a numerical approximation of A* (T) from 

A*{T) ~ A*(T) + (7a*(T). 

This means that we can approximate a path dependent option by using A* and 
an error term. 

Although this asymptotic approach is only valid for calculating the value of 
long maturity options, it is beneficial because 

1. We can get the value of A*(T) instantaneously since a* is written as 
elementary functions in most cases. 

2. We can apply this method to many path dependent options. 

3. We can apply this method to a large class of the underlying asset. 

As an example of the first case, when St is geometric Brownian motion, a* is 
given by and written as elementary functions. As an example of the second 
case, we give the result in the case of an Asian option, the Lookback option, and 
the barrier option. As an example of the third case, we give the result in case of 
a geometric Levy process, particularly for Brownian motion (BM) and for the 
Normal Inverse Gaussian process (NIG). In our method, the longer a maturity 
time is, the better the computation accuracy is. That is because our method is 
based on the asymptotic behavior of price as T tends to oo. This is a notable 
result since most methods, like Monte Carlo simulation, are effective when the 
variance of the underlying asset is small; i.e., the maturity time is short. 

This paper is organized as follows. In Section [21 we present a theorem which 
gives the principle of our method and give the short proof. In Section [3J we 
give some exmaples of a* when the underlying asset behaves as some popular 
Levy processes, and we see that these a* are written by elementary functions in 
these cases. In Section 21 we give two theorems which help us how to find the 
approximation function A* (T) when given A* (T). Also, we show that some path 
dependent options, such as an Asian option, a Lookback option and a barrier 
option, can be applied these theorems. In Section [SJ theorems of Section 2] are 
proved, and in Section [6J we present some properties of (BM) and (NIG) used 
in the proof. 



2 Generalized Principle 

We set the notations again. Let an underlying asset process, S s be a positive 
Markov process. For a fixed time r > 0, we denote the set of all cadlag functions 
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with domain [0, r] by D([0, r]). We also define an expectation of a path depen- 
dent option whose maturity time is T and whose monitoring period is [T — r, T] 



where g is a function from D([0, r]) into R. For given A*(T), we also define 



where g is also a function from Z)([0, r]) into KL In this paper, we treat this g 
as the approximation function of the path dependent option. For the sake of 
simplicity, we use the notation t := T — r, 



and we see the asymptotic behavior as t — > oo instead of the asymptotic behavior 
as T — > oo. Under these conditions, we derive the following theorem. 

Theorem 2.1. We assume two assumvtions, HOI and U^B about (g,g,S); 

Al. Process S decays uniformly; that is, there exist measures v and v and a 
funcion, a : [0, oo) — > (0, oo), such that 



by 



A*{T) := E[g((S s ) T -T<s<T)] 



A*(T) by 



A*{T) := E[g((S s ) T -T< s <T)\ 



A(t) := A*(T) = E[g((S s ) t < s < t+T )}, 
A(t) :=A*(T) = E[g((S s ) t < s < t+T )}, 




where v t {M) : 
continuous, an- 




A 2. 




where 



e(x) := E[g((S*)o< s < T ) - ff((^)o< s <r)] 



anrf S 1 ^ is a process starting at x; i.e., S\ 




Then it follows that when t — ► oo, 



(2.1) 
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We remark that, for a*(T) := a(t + t), since 



A*(T)-A*(T) _ A(t) - A(t) 
a*(T) a(t) 

CL2) means (|2~Tj) . 

Proof of the theorem. From the Markov property of S, we obtain 

A(i) - A(t) = E[g((S s ) t < s < t+T ) - g({S s ) t < s < t+T )] 

= E[E[g((S s ) t < s < t+T ) - g((S s )t< s < +T )\S t }} 
= E[e(St)). 

Since J \e(x)\ sup ( ^r{x) v{dx) < oo, the Lebesgue convergence theorem implies 
that 



a(t) J dv 



e(x) — (x) v(dx) 
e(x) v(dx). 

□ 



3 Asymptotic Order of each model 

In this section, we consider I All in Theorem 12.11 The existence of a(t) depends 
only on each process, S. It is clear that a is unique in the sense of order; i.e., 
a(t) ~ a'(t) as t — > oo if a and a' satisfy |AT1 for the same process. We give 
some examples of St satisfving lAll represented as St = e Zt , where Z t is a popular 
Levy process. In the following argument, we regard measure v as a Lebesgue 
measure. 

For the sake of simplicity, we discuss Z instead of S. Let vt be a modified 
distribution of Z t \ i.e., Vt{A) := „m ■ Then, since St = e Zt , we have that 

dv tl , d P(Z t < z) d P(e Zt < e z ) de z d P(S t < e z ) dv t , zx 

(z) = - - = - - = - - = e (e ). 

dv dz a{i) dz a(t) dz de y a(t) dv 

Thus, by substituting e y — x, the condition [XT] is replaced by a condition about 
Z; 

Lemma 3.1. v t ^S> v, v ^S> v, and 



for any z imply EH for %{x) := i^f(logx). 



5 



3.1 Brownian Model 



Let us define Z t = z + aW t + [it where Zq, /i £ K, a > 0, and {W t }o<t< o is a 
one-dimensional Brownian motion. Then 

, . 1 . n 2 t. 

a{t) = _ exp( __) 

and % {z) = v^ exp O 

satisfy (pHj) . 



3.2 Normal Inverse Gaussian Model 

Let us define Z t — z + Wjq,^ + 9IG(t) + bt where z Q , b e R, and we denote 
the inverse Gaussian subordinator, IG, by IG(i) = inf{s > 0; B s + (is > St}, 
where /i e R, { > 0, and B is another Brownian motion independent of W. 
This process is called normal inverse Gaussian, see |S]. Then 

a(t) = t"5 exp (t(uS -9b- ^(b 2 + 5 2 )(d 2 + m 2 ))) 



<^ v / 2tt(6 2 + S 2 ) \b 2 + S 2 



satisfy (pHj) . 



3.3 Variance Gamma Model 

Let us define Z t = z Q + aW^u\ + 0j(t) + mt, where z > 0, 0, m e R, and we 
denote the gamma process by j(t) with variance rate A, see |16j . Then 



f 1-3 8m w 

t — x 



, dv. , 1 / 2 + A0 2 /^ 2 A ?7 - 1 w , 

md Tv [z) = 2 V 2^(1+^ eX H ( ^ + ^T )(Z " Z0) 



satisfy (E3|) where tj := x /l + + £). 



4 Path dependent Payoff and its Approximation 
function 

In this section, we consider IA2I in Theorem 12.11 To satisfy lA"2l for (g,g,Z), we 
need to select a "good" g for a given path dependent payoff, g. The following 
theorems give classes of (g,g), which satisfy The first class includes, as 
examples, Asian options and a Lookback option. 
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Theorem 4.1. Let h : D([0,t]) -> R satisfy 

(4.1) inf w(s) < h(w) < sup w(s), 

0<S<T 0<S<T 

(4.2) h(aw) = ah(w), 

for any w £ -D([0, r]) a«c? a € R, and Zei (g, <f) represent 

g(w) = (h(w) - K)+, g(w) = C(e w ^ - K/C)+, 

where C :— E[h((St)o<t<T)]- Then \AB holds in cases o/(BM) and (NIG) with 
the parameter condition, 

— »_<!. 
(5 + V6 2 + 5 2 

The second class includes a barrier option. It also includes the Asian options 
with another g. 

Theorem 4.2. Let (g,g) satisfy there exist Mjj> Ml and M such that 

(4.3) inf w(s) > M L =>■ g(w) = g(w), 

0<s<r 

(4.4) sup w(s) < Mu g(w) = g(w), 

0<S<T 

(4.5) \g(w) - g(w)\ < sup w(s) + M. 

0<S<T 

Then VAM holds in cases of (BM) and (NIG) with the parameter condition, 

26 <1. 



S + v*b 2 + S 2 

We note that the classes given by Theorem 14.11 and Theorem 14.21 are just 
examples of class satisfying IA2I and that there may exist many other classes. 
Especially, for given payoff g, the g is not unique. 

4.1 Discrete Asian option 

The payoff for a discrete Asian option whose strike price is K and whose matu- 
rity time is T is defined by 



where < t\ < ■ ■ ■ < tjv = t. In this case, let h be defined by 

JV 

A' 



1 N 

K w ) ■= mz2 w ( t ~ t 0- 
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We can then express the option as 

g((S s ) T -T< s <T) = (h((S s ) T -T<s<T) - K) + 

Then it is easy to see that h satisfies the conditions of Theorem 14.1 
To apply Theorem 14. 2\ let (g,g) be defined by 

ffH := ( T?) ]w(T-Ti)-K\ , 



\ i=i 



i N 

i=l 

where Yli=i ^ = Then (g, g) satisfys the conditions of Theorem [521 with 
respect to Mjj := min^ Ki and Ml := max^ Ki. 

4.2 Integral Asian Option 

Similar to the payoff for a discrete Asian option, we define that for an Asian 
option as 



t > 
S s ds — K 

T-t j 



In this case, let h be defined by 



1 r 

h(w) .= — I w(s)ds. 
T Jo 

Then h satisfies the conditions of Theorem 14.11 
To apply Theorem 14. 2\ let (g,g) be defined by 

g(w) := ^— J w(s)ds — K 

g(w) ■= - [ M«) - k s ) + ds 

T Jo 

where \ J Q r K s ds — K. Then (g,g) satisfys the conditions of Theorem 14. 21 with 
respect to Mjj := inf s K s and Ml '■= sup s K s . 

4.3 Discrete Lookback Option 

The payoff for a discrete Lookback option whose strike price is K and whose 
maturity time is T is defined by 

max St—t- ~ A' 

=1,...,N 



where < t\ < ■ ■ ■ < tjv = t. In this case, let h defined by 

h(w) := max w(tA. 

y ' i=l,...,N ^ ' 

Then h satisfies the conditions of Theorem 14.11 
4.4 Partial Barrier Option 

The payoff for a knock-in partial barrier option with the strike price K, the 
barrier level L, the maturity time T, and monitoring period [T — r, T] is defined 

by 



(St-K^I inf S t >L} 

T-T<t<T 

where L > K. Let (g,g) be defined by 

g(w) := (w(t) - K) + l {infw > L} , 
g(w) := (w(t)-L)+. 

Then (g,g) satisfies the conditions of Theorem 14.21 with respect to Mjj := K 
and M L := L. 

5 Proof of Theorem 14.11 and Theorem 14.21 

To prove theorems, we use Lebesgue integrability. By substituting e z = x, the 
condition I A2I is replaced by a condition about Z; 



Lemma 5.1. 



/OO 
\e(z) | sup — (z) v{dz) < oo 
- oo t 



implies \AS\ for e(z) :— e(e z ). 

To see (|5.ip . it is enough to demonstrate that the integrand vanishes fast 
when \z\ — > oo. At first, we start with descriptions of the bound of sup t ^j-(x). 
Note that this density function, sup t ^-(x), depends only on the process Z but 
not payoff and its approximation function, (g,g); 

dv, 



(5.2) (BM) sup -£(z) < C ( BM) exp (JL Z ) 
t dv Vo" 2 / 



(5.3) (NIG) sup § (z) < C mG) exp ( 1 6 + 



2; 



t 



where C(bm) an d C(nig) are constants. 



9 



We next give a lemma about the order for 

e(z) = E[g{{e z °) Q < s < T ) -g{{e z ") < s < T )] 

to prove Theorem 14.11 where Z z s :— Z s — zq + z. The order is described by the 
tail order for each process, Z. 

Lemma 5.2. Let h and (g,g) satisfy the same conditions as in Theorem \4-l\ 
Let Z® be a Levy procces starting with such that 

< min{ inf P{Z° U > 0) inf P(Z° U < 0)}. 

0<u<r 0<u<r 

Then it holds that as z — » oo ? 

\e(z)\ = O (P(Z° T < log K-z)). 
Also, for any e > 0, it holds that as z — > — oo, 

\e(z)\ = 0(P(Z°> log K-zf- e ) 

ifE[ei supZ °] < oo. 

We next give a description of the order when Z is either (BM) or (NIG); 

lQgiC + M T|) 2 

2o- 2 x J [Z — > OOJ 

log Jf + HT|) 2 . 

2<7 2 T J — > — OO) 



(5.4) (BM) 


P(Z t < log K-z) 


= o(V i 

v z 


(5.5) (BM) 


P(Z t > log K -z) 


= oC-e- L 

z 


(5.6) (NIG) 


P(Z t < log K -z) 


= 0(\z\--* 


(5.7) (NIG) 


P(Z t > \ogK-z) 


= 0(\z\-i 


We now prove Theorem 14. 11 





) (z^-oo) 



Proof of Theorem \4-.l\ In the case of (BM), we have the following from the 
results of Lemma (|5lfl) . (O, and ([53)1 that 

|e(^)|sup^(z) = ofe~^ 

for any p < 1. This guarantees the integrability of the integrand in (|5.1[) . 

In the case of (NIG), we also have the following from the results of Lemma l5.2 
((Of . (j5T6|) . and that 



\e(z)\ sup — (z) = 0(|z| 2 ^e s+v^+a 2 v ) 

t <w 

for any p < 1. If 

26 

this guarantees the integrability of the integrand in (|5.ip . □ 
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Proof of Theorem Similarly to Thcorem l4.11 the theorem is proved by using 
a following lemma instead of Lemma 15.21 □ 



Lemma 5.3. Let {g,g) satisfy the same conditions as in Theorem \4-S\ Let Z® 
be a Levy procces starting with such that 

< mm{ inf P(Z° > 0), inf P(Z° < 0)}. 

0<U<T 0<1A<T 

Then for e > 0. it holds that as \z\ — > oo, 

\e(z)\ = O (mm{P(Z Q T > log M v - z) 1 ^ , P(Z° T < log M L - z) 1 ^}) 

l/jB [ e i(supZ»-infZi?)] < ^ 

5.1 Proof of lemmas 

Note that Lemma 15.21 and Lemma 15.31 do not require that Z be exact for either 
the (BM) or the (NIG) case, but that Z has a "reflection principle" ; 

Lemma 5.1. Let Z® be a Levy procces starting with 0. Then for any a > 0, 
P(Z° > a) > DP( sup Z s ° > a) and 

0<S<T 

70 ^ „\ ^ nut ;„t ?0 



wher 



P{Z»<-a) > DPI mi Z»<~a). 

0<s<r 



D := min{ inf P(Z° > 0) inf P(Z° < 0)}. 

0<U<T 0<U<T 



Proof. Let r a be the hitting time of (a, 00); i.e., r a := inf{s | G (a, 00)}. 
bmce the path of Z° is right continous, Z® > a if r a < 00. Therefore, the 
strong Markov property implies 

P{Z°>a) = P{Z°>a,T a < T ) 

> P(Z°>Z° Ta ,T a <T) 

= P(Z° T - Z° Ta > 0, r Q < r) 



[ T P(Z a T _ s >0)P(r a = ds) 
Jo 

inf P{Z Q U > 0) / P(r a = ds) 
' U < T Jo 



> 

0< 

= inf P(Zl > 0)P(r a < t) 

0<U<T 

= inf P{Zl > 0)P( sup Z° > a). 

0<U<T 0<S<T 

The proof of the other case is similar. □ 
We now prove Lemma 15.21 and Lemma 15.31 
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Proof of Lemma \5.2i From the definition of e(z) and the assumption of (g,g), 

e(z) = E[g((e z °) a < s < T )-g((e z *) < s < T )} 

= E[(h((e z "° ) < S < T ) - K)+ - C(e z « - K/C)+] 
= E[(h((e z+Z ") < S < T ) - K)+] - C(e z - K/C)+. 

When z < log ^ , we have the following inequality from (|4.ip , the Holder in- 
equality and the reflection principle that 

e(z) = E[(h((e z + z °)o< s <r)-K)+] 

< E[{ sup e z+z " - K)+] 

0<S<T 

< £[(e 2 exp( sup Z a s ) - K) p ] 1/p P{ sup Z° s > \ogK - z) 1/q 

0<s<t 0<s<r 

< E[(K/Cexp( sup Z Q S ) - /f) p ] 1 / p P( sup Z° s > \ogK - z) 1/q 

0<S<T 0<S<T 

< D-^ q E[{K/Cexp( sup Z° s ) - K) P ]^ P P(Z° > logK - z) 1/q , 

0<S<T 

where 1 < q < oo with i + A = 1. This yields the result 
e(z) = O (P(Z° >\ogK- z) 1 '^ 

as z — > — oo. 

On the other hand, when z > log |y, (|4.2p . (|4.1[) and the reflection principle 
imply that 

e(z) = El(h((e z *)o< s <r)-K)+]-(Ce z - K) 

= E[(h((e z+Z °) Q < S < T ) - K)+] - (E[h((e z ") Q < s < T )]e z - K) 
= E[(e*h((e z ») < s < T ) - K)~] 

< E\(e z inf e z "-K)-} 

0<S<T 

< KP{ inf Z° <\ogK - z) 

0<s<r 

< ^P(Z° T < log K-z). 

This yields the result i(z) = 0(P(Z® < \ogK — z)) as z — > oo. This completes 
the proof. □ 



Proof of Lemma \5.3[ Note that 



' - e*irde z ^— infe z e^ Z ^ z . 

inf e 



sup e s = e sup e - — ^ o - - z0 
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= E\ 



< E[(supe z * +M) 1, 



Then from the definition of e(z), (|4.3|) . (14.4)) and (|4.5j) . we have 
e(z) = S[.g((e z ?)o< s < T )-5((e^)o< s <r)] 

; [(3(( eZ Oo< S <r) - g((e Zs )o< S <r)) l{supe z J>M [/ ,M i >infe z f}] 
^{supe z ?>M[/,Af L >inf e z i }\ 
E\(M L ^ Z °-^ Z °+M) l {supeZ f >Ma , Mi >infe-f } ]- 

Therefore, the Holder inequality and the reflection principle imply 

e(z) < (M L E[e p{supZ °- [nl z '^]p + A/)P(supe z = > M v , M L > inf e z " )i 

= (M L P[e^ sup z "- inf z ")] p + M )P(sup Z° > log M v - z, log M L - z > inf Z°) * 

< (M L £:[e p(supZ °- inf z = 0) ]p + M)D _ i min{P(Z° > logAf^ - z), P(logM L - z > Z°)}i 

□ 

6 Appendix 

In this section, we present some properties of (BM) and (NIG); a and v of 
each process, ([572")l . (|573")l . (154]) . ([575)1 . ([576]) and (577)1 . 



6.1 Brownian Motion 

Recall that Zf = z + aW t + [it and Z t = Z%° . Then the density function of Zf 
with respect to the Lebesgue measure is 



y 1 



1 (a-pt-z) 2 

2cr 2 t 



\f2-naH 



By setting a : [0, oo) — > [0, oo) as a(i) = t - 2 e 2^ , we have 
— -(z = -^=e ^^+— ?2— . 

From the fact that < e 2,7* t < j f or a n ^ we h ave 

(ij>t . . 1 y('-'o) 

sup— — (z) < e ^ . 



This inequality shows (|5.2|) 
For y > 0, we have 



P(|i*|>|,) = *(^) + *(^), 



13 



where 

f°° 1 *2 

It holds from L'Hopital's theorem that 1 - $(£) = $(£) ~ -^e"V. This 
shows (El) and iPi 



6.2 Normal Inverse Gaussian 

Recall that Z* = z + WjQ, t s + 8lG(t) + bt where fi G M, and we denote the 
inverse Gaussian subordinator IG(i) = inf{s > 0; B s + fis > St} where S > 0. 
Then the density function of Z^ with respect to the Lebesgue measure is 



1 2 +^ 2 g/rft+flfc-.-w)^ ^ / (02 + M 2 )(1 + f 

*^ +1 V V V 7 

where K\ is the Bessel function of the third kind (See [5]). 

We use a estimation of the Bessel function from page378 of pQ; 



2 



(6.1) KM ~ y-e 

We set as 

a ( f ) = J_ e t(M-5-e6-V(6 2 +<5 2 )(e 2 +M 2 )) 



Then it holds that 



di>t , . ^ 1 / 6> 2 + /i 2 
sup— (z) < 



dv ~ J2n V ^ 2 



V fe+fc 2 v^+^ V^ Zo) 
e V a+%/i> 2 +s 2 ) 



This shows (IST2 



Note that St\j (9 2 + ^ 2 )(1 + (^r^f) -> oo when |z| ->• oo. By using 
6.1|) again, we obtain that the density function of Z® decays in the order as 



when \z\ — * oo. By using L'Hopital's theorem, we have (|5.6p and (J57TJ). 
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